Adsorption on carbon nanotubes: quantum spin tubes, magnetization plateaus, and 

conformal symmetry 
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We formulate the problem of adsorption onto the surface of a carbon nanotube as a lattice gas on a 
triangular lattice wrapped around a cylinder. This model is equivalent to an XXZ Heisenberg quan- 
tum spin tube. The geometric frustration due to wrapping leads generically to four magnetization 
plateaus, in contrast to the two on a flat graphite sheet. We obtain analytical and numerical results 
for the magnetizations and transition fields for armchair, zig-zag and chiral nanotubes. The zig-zags 
are exceptional in that one of the plateaus has extensive zero temperature entropy in the classical 
limit. Quantum effects lift up the degeneracy, leaving gapless excitations which are described by a 
c = 1 conformal field theory with compactification radius quantized by the tube circumference. 

PACS: 67.70+n, 7f.lO.Pm, 75.fO.Jm, 



Monolayer adsorption of noble gases onto graphite 
sheets has proven to be an interesting problem both theo- 
retically and experimentally . Many of the observed 
features can be understood within a lattice gas model, 
where the underlying hexagonal substrate layer forms a 
triangular lattice of preferred adsorption sites. An equiv- 
alent formulation is in the language of spin models on a 
triangular lattice, where the repulsion between adsorbed 
atoms in neighboring sites translates into an antiferro- 
magnetic Ising coupling. The frustration of the couplings 
by the triangular lattice leads to the rich phase diagram 
of the monolayer adsorption problem Q. Introducing 
hopping adds quantum fluctuations, further enriching the 
phase diagram 

In this paper we address what happens if, in addition 
to the triangular lattice frustration, one has an extra ge- 
ometric frustration due to periodic boundary conditions. 
In fact, such a system is physically realized by a single 
walled carbon nanotube which may be viewed as a 
rolled graphite sheet. In this context, adsorption has 
been the subject of growing experimental and theoreti- 
cal interest spurred by potential applications. Stan 
and Cole have considered the limit of non-interacting 
adatoms at low density, finding that they are localized ra- 
dially near a nanotube's surface at a distance comparable 
to that in flat graphite (~ sA). In that work, it was suf- 
ficient to omit the hexagonal structure of the substrate. 
However, the corrugation potential selects the hexagon 
centers as additional commensurate localization points 
0. In view of the similarity to fiat graphite, we include 
both the substrate lattice and adatom interactions and 
consider a wider range of densities. In fact, very recently, 
it has been shown |^ that the adsorbate stays within a 
cylindrical shell for fillings less than « 0.1/A (or « 0.5 
adatom/hexagon), justifying the densities studied here. 

The adsorption sites thus form a triangular lattice 
wrapped around a cylinder. The geometrically inequiv- 
alent ways in which the wrapping is realized are labeled 



by two integers, {N, M); the case {N, N) is known as the 
armchair tube, {N, 0) is the zig-zag, and all others are 
chiral. The geometric frustration is present whenever the 
wrapping destroys the tripartite nature of the triangular 
lattice, which occurs when {N — M)mod 3 is non-zero 
(this criterion is familiar in the context of electronic con- 
ductivity 1^). Thus, both zig-zag and chiral tubes can 
be frustrated geometrically, whereas armchair tubes can- 
not. In Fig. |l| we show the tube obtained from the (7, 0) 
zig-zag. In terms of the equivalent spin systems, this is 
similar to recent models of "spin tubes" . 




FIG. f. Adsorption sites on a (7,0) zig-zag nanotube 

When the adsorbed gas is a hard-core boson, the lattice 
gas is defined by the Bose-Hubbard Hamiltonian ||3|,|l0|| 
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where Ui is the boson density at site i, V is the nearest 
neighbor repulsion and t is the hopping amplitude. In 
the equivalent Heisenberg spin representation, 

n=-2tY,s^s^+s!sy + vY,sts;r-HY,s! , (2) 

where S'f = ni — 1/2 and H — fx — is an effective 
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external magnetic field. Throughout the paper we will 
use the spin and density representations interchangeably. 

The Ising limit, in which hopping is not allowed, al- 
ready contains many interesting features. We start the 
analysis in this regime, obtaining the phase diagram as 
a function of the magnetic field, and then consider quan- 
tum fluctuations perturbatively in t/V. We summarize 
our results first. 

The phase diagram in the temperature-magnetic field 
plane of a typical tube is shown in Fig. 0. 
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FIG. 2. Phase diagram of the (iV, M) tube 

When the index q — [N — Af )mod 3 is 1 or 2, we find 
four lobes (solid lines), corresponding to two plateaus 
with magnetizations m_ < 1/3 and m+ > 1/3. Here, we 
use the standard Ising notation in which spin is ±1. Note 
that the plateaus are real phases only at zero temperature 
because the tube is one-dimensional. At finite tempera- 
ture, the boundaries should be interpreted as crossovers. 
Nonetheless, deep within a lobe, at fc^T ^ V , the mag- 
netizations are well-defined. Specifically, for g = 1, we 
obtain the exact expressions 
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The complementary case of q = 2 is obtained by in- 
terchanging N ^ M. On the other hand, those tubes 
without geometric frustration {q = 0) behave similarly 
to the fiat sheet (dotted lines) which has only two lobes 
with magnetizations ±1/3 [||. As the tube perimeter 
approaches the flat sheet limit, one expects that the geo- 
metric frustration becomes irrelevant. Indeed, as N or 
M —>■ GO, m+ and m_ squeeze 1/3 as the inverse of 
the tube diameter and become indistinguishable. Beyond 
the lobes, where the field is strong enough to overcome 
all nearest neighbor bonds {\H\/V > 6 at ksT — 0), 
the tube is fully polarized. The filling fractions are ob- 
tained from the magnetizations by m — — 2(n — 1/2). 
The phase diagram, however, is more easily visualized in 
terms of spin since spin reversal, m ^ —m, corresponds 
to particle-hole symmetry, n ^ 1 — n. 

We have verified this prediction numerically by transfer 
matrix methods for zig-zag tubes up to = 11 and for 
the chiral tubes up to N+M = 7. Although 7 is probably 
too small to be physical, we believe that the arguments 



in this paper generalize to any tube. In Fig. |^ we display 
sample data for two zig-zag tubes with different q: (7, 0) 
and (8, 0). The magnetization curves show clear plateaus 
whose values and transition fields match those predicted 
by Eq. (||). By increasing the temperature and follow- 
ing the evolution of the plateaus, we generate the phase 
diagram above. 

We find that a rather interesting feature of the zig-zag 
{N, 0) tubes emerges, making them exceptional. The in- 
sets in Fig. ^ indicate an extensive entropy at zero tem- 
perature, which has plateaus, too. Upon enumerating 
the degenerate space explicitly, we shall show that the 
entropy is exactly s = (In2)/A^ and that it occurs in m+ 
for q = 1 and in m_ for q = 2. In the presence of 
hopping, the non-degenerate plateaus retain their gaps, 
whereas the degenerate ones become correlated states 
with a unique ground state and gapless excitations. More 
precisely, conformal invariance develops and the effective 
theory has central charge c = 1 with a compactification 
radius, R, quantized by the tube circumference, R = N. 
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FIG. 3. Magnetization and entropy per site at ksT — 0.05V 

In order to understand the magnetizations and nature 
of the geometric frustration, it is more intuitive to use 
the original bosonic picture. As a result of hard-core 
repulsion on the infinite graphite sheet, the m = 1/3 
plateau corresponds to filling one of the three sublattices, 
A, B or C, of the triangular lattice. This configuration 
minimizes the repulsion, VniUj, while maximizing the 
filling, /in. It is natural to try the same for nanotubes, 
as we illustrate in Fig. ^ for (5, 0). 




FIG. 4. Fillings and zipper of the (5, 0) zig-zag tube. n± 
corresponds to m^i 
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Upon wrapping, however, the thick vertical hnes are 
identified and the lattice is no longer tripartite. In fact, 
the number of sublattice sites is no longer equal, and 
there is a mismatch along the thick line, which we term 
the "zipper" . On the left we fill the A sublattice, obtain- 
ing the filling fraction ri+ = 2/5, and on the right either 
B or C may be filled with the result that n_ = 3/10. 
For general (A^, 0) there are 2N hexagons in the unit 
cell, and the filling fractions are n+ = \2N /'S\/2N and 
ri_ = [2A^/3J/2A^, where \x~\ and [a:J denote the larger 
and smaller of the two bounding integers of x, respec- 
tively. The magnetizations in Eq. (^) follow directly by 
using the correspondence m = — 2(n— 1/2). Furthermore, 
due to the sublattice mismatch, the number density of 
adjacent particles, rib, may be non-zero. In the case of 
(5, 0), there are two broken bonds per unit cell in n+, and 
none in 7i_. This result generalizes to any q = 2 zig-zag 
tube: = 2/2iV and = 0. For q = 1, the argu- 
ment goes through as before, except that = 1/2N. 
We summarize this compactly by nf,+ — q/2N. 

Substituting these fillings into the Hamiltonian (|l|) 
yields two energies per site, e±{^) = Vni,±— ^n±. The 
transition occurs when these levels cross: e+ = e^, or 
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Solving for /i and using the correspondence H = fi — 3V 
gives precisely the critical field in Eq. (^). In particu- 
lar, this explains why there are exactly two independent 
plateaus. Note that, for the special case of the zig-zags, 
the critical field depends only on q and not on N per se. 

In the above analysis, we have made only one assump- 
tion, namely that the zipper runs parallel to the tube 
axis. In general, the zipper may wind helically around 
the tube or wiggle sideways. However, in all the cases 
that we considered, the straight zipper has the lowest 
energy, and moreover, our transfer matrix computations, 
which are blind to this assumption, arc consistent with 
our analysis. 

The chiral tubes are different. Due to their geometry 
the zipper is forced to wind, but, again, we find that the 
choice of the straightest possible zipper reproduces our 
numerics for N + M up to 7. The determination of the 
fillings and level crossings is much more involved than 
that of the zig-zag, and we leave it for a more detailed 
paper. In any case, our analysis reveals that the plateaus 
in a chiral tube are not macroscopically degenerate, so 
that the zig-zags are at a special degenerate point. 

Having understood in detail the Ising limit, we now 
turn on a small hopping, t V ^ that introduces quan- 
tum fluctuations. Deep within a plateau, the substrate 
is maximally filled since adding a particle increases rif,. 
Consequently, all plateaus begin with a classical gap of 
order V , and we work in the Hilbert space of the classical 
ground states. Those plateaus which have only a discrete 




symmetry must retain their gaps, but the macroscopi- 
cally degenerate plateaus are more complicated. 

Let us reconsider the n+ fiUing of the (5,0) tube in 
Fig. ^. Notice that a particle may hop laterally by one 
site without changing n;,, as we illustrate in Fig. ||, left. 
Imagine building a typical state layer-by-layer from 
top to bottom, with a total of L layers. Each new layer 
must add exactly two filled sites and one nearest-neighbor 
bond (rif, — 1/5). This constraint implies that no two 
adjacent sites may be occupied within a layer; if they 
were, then, to conserve rib, two adjacent sites must be 
occupied in the next, and so on up the tube. However, 
this state is not connected to any other by a single hop. 
Similarly, the particles cannot hop from layer to layer 
because this adds another intra-layer bond. 
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FIG. 5. LEFT: Typical configuration in n+ (or m_) of 
the (5, 0) tube. Alternating numbering within layers allows a 
symmetric description from bottom-to-top or top-to-bottom. 
RIGHT: Allowed states as paths on a wrapped square lattice. 
The vertex labels may be dropped. 

An allowed state can be represented as a string of oc- 
cupied sites, {cTi}, i = 1, . . . , L, which in our example is 
{•■•(5,3)(5,2)(3,1)(2,4)(1,4)(1,4)...}. At each layer, 
there are exactly two possibilities for the following one. 
For example, (1, 4) can be followed by (1, 4) or by (2, 4). 
However, the total number of possibilities at any given 
level is five. Fig. ^ (right) summarizes this structure suc- 
cinctly as a square lattice wrapped on the cylinder. A 
typical state, then, is a lattice path along the tube. 

Generalizing to (A^, 0), we find N possible states in 
each layer and two in the succeeding one, and the struc- 
ture of states is again that of a wrapped square lattice 
with N squares along the circumference. The dimension 
of the Hilbert space is the number of lattice paths, A^2-^, 
so that in an infinitely long tube, the entropy per site 
is exactly (ln2)/iV, as claimed earlier. Notice that con- 
strained paths introduce correlations along the length of 
the tube, despite the absence of inter-layer hopping. 

The matrix elements of the projected Hamiltonian con- 
nect only those states that differ by a single hop: 



(MI^IW) = { 



-2t if Y.^ 5a^r^ 

otherwise 



(5) 



We diagonalize this Hamiltonian numerically with peri- 
odic boundary conditions for system sizes up to = 11 
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and L = 10. Additionally, we can obtain the ground 
state energy up to L = 16 due to the sparseness of H.. 
We will fix 2t = 1 in what follows. 

We find that the degeneracy is lifted and the ground 
state becomes unique and uniform. The ground state en- 
ergy, Eo{L), follows Eo^-QmiL-Tic/Q with c- 1.005. 
The lowest — 1 excited states are given by = 
a'^A/{N'^L), with A = 12.9±0.5, which is shown in Fig. | 
for a — 1,2, 3. All of these levels are doubly degenerate. 




l/N^L (x 1000) 
FIG. 6. The gap scales as 1/N^L. 

This ground state energy and spectrum are in per- 
fect agreement with a conformally invariant bosonic 
theory with central charge c — I compactified on a 
radius R = C^- We take the Lagrangian density 
C — -^[v^-^ (dtcf))'^ — v{dx<p)'^], where u is a velocity, and 



the compactification is defined by = 
mode energies of this theory are ||ll[ 



E 
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2TrR. The zero 
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where a, b are integers that label topological momenta 
and windings of 4>. Right- and left-moving oscillator 
modes of energy ujn = vkn, where fc„ = 2Trn/L, also 
appear in the spectrum, but for a < N the zero modes 
are the lowest. Our spectrum in Fig. |6| corresponds to 
E^ J, with b — 0. To fix (, we look at higher low-lying 
levels (which also scale like 1/L). We find that the A^'th 
excitation energy is independent of A^ and quadruply de- 
generate. This can happen only if the A^'th zero mode, 
^±N ~ 27ru/C^L, is degenerate with the lowest oscilla- 
tor mode, U!±i = 2t:v/L, which fixes C 1- Thus, the 
compactification radius is R = N. The velocity can be 
read off from the slopes in Fig. ^ as u = A/27r. The rest 
of our spectrum is consistent with these parameters. 

One observable consequence of conformal symmetry is 
that the low temperature heat capacity is fixed by c 



C = 



-T 



•y I S 
0: 



(7) 



It is noteworthy that, even though the dispersion of 
the oscillator modes is independent of A^, the spectrum 
remembers, via the zero-modes, the finite radius of the 



nanotube. Furthermore, R is quantized by A"; in the 
language of Luttinger liquids, this means that the Lut- 
tinger parameter is fixed by topology, similarly to the 
case of edge states in a fractional quantum Hall fluid 
p2[ , and in contrast to quantum wires (where the Lut- 
tinger parameter can vary continuously). Because there 
is no inter-layer hopping, (j) is tied to transverse, rather 
than to longitudinal, density fluctuations along the tube. 
We will present a detailed analytical derivation of the 
effective theory from the lattice in Fig. ^ elsewhere. 

Let us briefly view the spin tube as a quantum spin 
ladder to see if it yields the zero gap. A standard ap- 
proach is to use a Lieb-Schultz-Mattis (LSM) argument, 
in which the spins are deformed slowly along the length 
|l3| . Applying it to our tube, we find that a plateau is 
gapless if S* — M is not an integer, where 5* is the total 
spin and M the magnetization per layer. Using S = N/2 
and the magnetizations from Eqn. (||), we find that S — M 
is an integer in the macroscopically degenerate plateaus, 
so that the LSM argument is insufficient in this case. A 
conclusive argument must take the geometric frustration 
into account. 

Before concluding, we should point out that the geom- 
etry of the (2, 0) tube is special; all sites in adjacent layers 
are interconnected. As a result, all of its plateaus have 
an extensive entropy, and we find that hopping opens a 
gap in both plateaus. In fact, this tube can be written as 
a spin chain that has been studied at isotropic coupling 
[Q, —2t = V. Two plateaus were found in this case, and 
it is tempting to speculate whether the two regimes are 
connected adiabatically. 

In conclusion, we have studied the problem of mono- 
layer adsorption on carbon nanotubes and identified sev- 
eral interesting filling fraction plateaus. Since the differ- 
ence between the plateaus decreases slowly, as the inverse 
of the tube diameter, experimental measurement should 
be feasible for large enough tubes. We have identified 
the zig-zag tubes as exceptional, in which the geometric 
frustration together with quantum fluctuations lead to 
conformal symmetry. This system is a physical realiza- 
tion of quantum spin tubes. 
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